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Abstract. We compute the anomalous dimension of the fermion field with flavours in the 
fundamental representation of a general Lie colour group in the non-abelian Thirring model at 
■ four loops. The implications on the renormalization of the two point Green's function through 



the loss of multiplicative renormalizability of the model in dimensional regularization due to 
the appearance of evanescent four fermi operators are considered at length. We observe the 
appearance of one new colour group Casimir, d'^'^'^d'^'^'^, in the final four loop result and discuss 
its consequences for the relation of the Knizhnik-Zamolodchikov critical exponents in the Wess 
Zumino Witten Novikov model to the non-abelian Thirring model. Renormalization scheme 
changes are also considered to ensure that the underlying Fierz symmetry broken by dimensional 
regularization is restored. 
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1 Introduction. 



The non-abelian Thirring model, (NATM) , is a two dimensional renormalizable fermionic quan- 
tum field theory with a variety of interesting properties, For instance, it is asymptotically 
free, and being two dimensional its bosonized version is related to the Wess Zumino Witten 
Novikov, (WZWN), model which is a bosonic nonlinear a model with a topological term or 
equivalently a torsion potential, Indeed its connection with the WZWN model has played a 
fundamental role in understanding two dimensional conformal field theories following the early 
study by Dashen and Frishman, of examining the conformal properties of the NATM and 
its fixed point structure. An interesting and elegant feature to emerge from the conformal ap- 
proach to solving the NATM or equivalently the WZWN model was in the work of Knizhnik 
and Zamolodchikov, Q , where the critical exponents of the fields and parameters of the theory 
were written down exactly. These all orders results depended in a simple fashion on the group 
Casimirs of the underlying non-abelian symmetry of the model. In particular the elementary 
Casimirs T(R), C2{R) and C2{G) were involved as well as Nj the number of fiavours if the fields 
had an additional internal symmetry. These exact results of Knizhnik and Zamolodchikov have 
been checked in a variety of ways. For instance, in |Q, ^ the perturbative renormalization group 
functions were computed in the nonlinear a model on the group manifold with a Wess Zumino 
term to several orders in perturbation theory. Then evaluating the expressions at the non-trivial 
fixed point the results were compared with the Knizhnik-Zamolodchikov exponents when these 
were expanded to the equivalent order in perturbation theory. In expressing the essence of such 
a check in this condensed way it is important to recognise the technical difficulties in perform- 
ing the renormalization of a theory with a topological term. For example, in [Q, |5[ the use of 
dimensional regularization led to the problem of loss of multiplicative renormalizability due to 
the appearance of an evanescent (kinetic) operator. Such operators only exist in d-dimensions 
and when restricted to two dimensions collapse or evaporate to zero. However, their presence in 
a dimensionally regularized renormalization is fundamental to ensuring the model can be ren- 
dered finite and they affect the determination of the true renormalization group functions in a 
subtle way. Ignoring their presence and effect would mean that agreement between perturbative 
expressions and the Knizhnik-Zamolodchikov results would not be possible as emphasised in 

ii- 

Having recalled these connections between the WZWN model and general results one nat- 
ural question which arises is to what extent can one derive or verify similar results in explicit 
perturbation theory in the NATM which is connected via the bosonization rules of |^]. To ap- 
preciate the non-trivial nature of such a problem it is worth considering a few salient properties 
of the fermionic model. One interesting feature which it possesses is a connection with quantum 
chromodynamics, (QCD), in d-dimensions. In Q it was argued that the strong coupling limit 
of the NATM in d-dimensions possesses a gauge symmetry. Moreover, if one examined the re- 
sulting effective action in d-dimensions in the large Nj limit where Nf is the number of quark 
fields, then the Feynman rules of QCD emerge naturally. In essence the closed quark loops with 
the requisite number of external points reproduced the Feynman rules for the triple and quartic 
gluon vertices in the limit of going to four dimensions, |Q. Higher point functions vanished in 
the same limit. This remarkable connection in d-dimensions with the NATM was exploited in 
a variety of papers, [^, to compute d-dimensional critical exponents. These provided all 

orders results in QCD perturbation theory in the strong coupling constant and agreed with the 
corresponding explicit high order perturbative results in the region of overlap. Moreover, one 
could go beyond the leading order in the large expansion and compute exponents at 0{1/Nj), 
1^], which agreed with the known four loop MS results of jl^, |ll|, so that the connection of |^] 
is confirmed to be not solely leading order. Since such a QCD NATM connection is now es- 
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tablished one can place properties of QCD perturbation theory in the context of the NATM as 
well as the WZWN model. For instance, the four loop /3-function of QCD has been computed 
in MS in |12|. One novel feature of this result was the appearance of new colour group Casimirs 
beyond the elementary ones already mentioned. They also arise in the quark mass anomalous 
dimension, |11], and must be present in the quark anomalous dimension at four loops. However, 
whilst this latter quantity has been calculated in QCD in |13|, it was only for the colour group 
SU{Nc). One cannot explicitly deduce the result for an arbitrary Lie group from this result since 
there is no unique construction of the new general Casimirs from the Nc dependent information. 
However, by examining the topology of Feynman diagrams it is clear that the new Casimirs 
can be present. Therefore, having recalled these connections and group theory one issue which 
arises is to do with the structure of the four loop anomalous dimension in the NATM. On the 
one hand the connection with QCD and the Feynman diagram topologies indicates that the new 
Casimirs ought to be present. However, the relation of the WZWN model and the Knizhnik- 
Zamolodchikov critical exponents suggests that when the renormalization group functions are 
evaluated at criticality either the new Casimir structures somehow cancel or they are not present 
in the first instance. Therefore, this is the main problem we address in this paper where we 
will compute the four loop fermion anomalous dimension in the NATM to ascertain the group 
structures which will appear. This extends the earlier two and three loop calculations of re- 
spectively 1 14] and |15]. However, it is far from a straightforward renormalization since we will 
use dimensional regularization and like the treatment of the WZWN model the NATM ceases 
being multiplicatively renormalizable. This was recognised in [16, 17, 18, 1£] and the necessary 
formalism was introduced to handle the contributions from the new evanescent four fermi op- 
erators which result in the renormalization based on pO|. In |l5| the dimensionally regularized 



two loop /3-function calculation of |14| was extended to three loops and the terms beyond that 
originally computed in the cutoff regularized calculation of Destri, were produced. Whilst 
the result for a general colour group was quoted for the three loop /^-function that calculation 
exploited the fact that to this loop order one could in fact compute in the case when the colour 
group was SU{Nc) and still construct the general result uniquely at the end. This was done 
to speed intermediate aspects of the computation. As already indicated since such a procedure 
would fail to address the Casimir problem at four loops, we choose here to consider a general 
(classical) Lie group and not specify to an SU (A^) group in contrast to |l|]. Whilst this will lead 
to a longer calculation it will reveal a rich group structure at various stages. Moreover, we will 
address the evanescent operator issue in the general group case and establish the renormaliza- 
tion group functions akin to those which underpinned the check of the Knizhnik-Zamolodchikov 
exponents in the WZWN model. Given the nature of the NATM being a four fermi interaction 
our calculation is intimately connected with the 0{N) Gross Neveu model, [21|, whose four 
loop anomalous dimension is known in MS, |^2|, which extended the lower order calculations 
of [p3| , p^ . This will provide an important cross-check on our integration routines which have 
been implemented in a symbolic manipulation package called Form, |25|. By ensuring that 
the correct four loop Gross Neveu result is established in our programmes, changing the input 
Lagrangian to the NATM means one can essentially rely on that aspect of our calculation being 
correct, aside from the usual internal checks one usually has. Although we will produce the MS 
result it is important to note that in the NATM, the Fierz symmetry of two dimensions gets 
destroyed when one computes with a dimensional regularization, |1(:, 17, 15, 19|. Therefore, we 
will need to consider a finite renormalization beyond minimal subtraction in order that the Fierz 
symmetry and other model equivalences are restored in the perturbative renormalization group 
functions similar to those considered in [^. For example, the single flavour abelian Thirring 
model is equivalent to the single flavour Gross Neveu model. However, by examining the three 
loop MS NATM fermion anomalous dimension, [^], the restriction of that result to compare 
with the same MS result in the Gross Neveu model shows that it is not in agreement and there- 
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fore in need of a finite renormalization to restore its equivalence. Whilst this additional finite 
renormalization may at first sight appear to be technical, it is in fact a standard way of dealing 
with the limitations of dimensional regularization. In other words it does not preserve certain 
symmetries of the theory, such as chiral symmetry, in four dimensional gauge theories. The 
main benefits of using a dimensional regularization, which far outweigh other approaches, is the 
ability to calculate to high loop order. 

The paper is organised as follows. In section 2 we recall the basic properties of the NATM 
necessary for the four loop renormalization including the relevant lower order renormalization 
constants. Those for the four fermi evanescent operators have been rederived for the case 
of a general colour group. Section 3 is devoted to the technical details of the full four loop 
renormalization including how certain two loop subgraph Feynman integrals are evaluated by 
the Gram determinant method of and how the evanescent operators are handled in the 
higher order corrections. The results of the full calculation are discussed in section 4 where a 
study of one general set of scheme changes that can be performed is given. Our conclusions are 
given in section 5 and an appendix is devoted to the derivation and general aspects of the Gram 
determinant method. 



2 Background. 

The Lagrangian for the strictly two dimensional (massless) non-abelian Thirring model is 

2^natm ^ + | (^^^-f ^M^^^^i^) ^ (2.1) 

where -0*^ is a Dirac fermion with flavour and colour indices i and / respectively with 1 < i < Nj 
and 1 < I < iVfyjjfj. Here we denote by A^f^j^^j the dimension of the fundamental representation 
of the colour group G whose group generators are T", 1 < a < N^, and is the dimension of 
the adjoint representation. The generators obey the usual Lie algebra 

[T^T*] = if^^bcrj.c ^2.2) 

where /"'"^ are the structure constants. We have chosen the sign of the coupling constant g in a 
non-standard way. This is in order to make extensive use of the results of |15| from the point of 



view of renormalization constants but it is elementary to set g = — X where A is the conventional 
coupling constant at the end of the calculations to ensure the model is asymptotically free, M. 



Unlike |15], however, we will not specify the colour group to be SU{Nc) since we will be concerned 
with determining the rich Casimir structure of the anomalous dimension at four loops. The 
advantage of restricting attention to SU{Nc) in determining the three loop /3-function was that 
the final result could only depend on the elementary group Casimirs C2{R) and C2{G) as well 
as T{R)Nf where T'^T'* = C2{R), f'"^/'""^ = C2{G)5''^ and Tr(T''T'') = T{R)6''K Therefore, by 
working with T{R) = i, Csfi?) = {N^ - l)/{2Nc) and C2{G) = for SU{Nc) it was possible 
to uniquely determine the result for general G from the SU{Nc) value of the /^-function. The 
potential appearance of new Casimirs at four loops means this avenue is closed to us here. 
Therefore, we are forced to work with a general Lie group. For definiteness the (three) new 



Casimirs are products of the symmetric tensors, [12 



where (^'^)bc = ~ if"-^'^ is the adjoint representation of the Lie algebra. In particular for SU (Nc), 

an, 

dfcd^abcd ^ (jv,4 _ 6Ar2 ^ 18) df^'^df^'^ _ Nc{N^ + 6) 

Aa ~ 96A^2 ' ]va ~ 48 ' ^ ' ^ 
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Whilst we have to adopt a different strategy to deal with the group theory of (2J) the 



approach to compute the anomalous dimension of remains the same as [|^, |l^, 17, 18 
In 1 14] the two loop calculation was performed in strictly two dimensions using a momentum 



cutoff and a version of (^2) which involves an auxiliary spin-1 field, which in d-dimensions would 
correspond to the gluon of QCD at the large non-trivial fixed point, Q. Indeed this is the 
reason for choosing ij)^^ to be in the fundamental representation of G. However, using a cutoff to 
regularize the integrals is not a feasible approach beyond two loops despite the advantages that 
exist from maintaining a 7-algebra which is strictly two dimensional. Instead it is better to use 
dimensional regularization retaining the MS scheme since it allows one to compute the massless 
Feynman integrals more easily. The major difficulty in this is the loss of the finite dimensional 
7-algebra. See, for example, ||2^. In two dimensions, for instance, the identity 

t'^t'Tm = (2.5) 

reduces the structure of the numerator of an integral substantially given the nature of the 
interaction of (^]^). In d-dimensions one must use 

l^Yl^. = - (d-2)7" (2.6) 

so that extra terms will emerge. However, there is a deeper subtlety involved and that is that two 
dimensional four fermi theories cease to be multiplicatively renormalizable in d-dimensions. For 



the NATM the breakdown is at one loop and was pointed out in [16, |l^. Latterly for the Gross 



Neveu model it has been shown in |29] that this model is also not multiplicatively renormalizable 



in dimensional regularization with the first occurence being at three loops. This has been 
verified in |15|. This feature is evidenced in the appearance of evanescent operators in the 



renormalization of the four point function. In other words new four point operators are generated 
under the renormalization which do not have the structure of the original interaction term and 
have the evanescent property that they do not exist in the two dimensional limit. This may 
suggest that these extra operators play no role in the construction of the strictly two dimensional 
renormalization group functions. However, if we focus on the NATM beyond one loop one has to 
include these new operators in the higher loop Green's functions and more importantly account 
for their effect and contribution in the renormalization group functions in d-dimensions prior to 
restricting the MS results to two dimensions. We recall that when one calculates the /3-function 
in a theory (without evanescent operators) one works with renormalization constants, Z, which 
depend on e where d = 2 — e and a coupling constant which is dimensionless in d-dimensions. 
In deducing f3{g) correctly the dimensionality of the coupling constant in d-dimensions is crucial 
in restoring terms which are of the form e/e. Therefore when evanescent operators are present 
they in fact subtlely contribute to the renormalization group function in such a way that naively 
following the usual procedure described above, their presence remains in such renormalization 
group functions. Hence one would obtain incorrect results for the strictly two dimensional 
functions. Indeed this point can best be appreciated if one calculates the two loop /3-function of 
the NATM using the naive approach. In a result emerges for this /5-function which does not 
agree with the two loop result of It is obvious that the d-dimensional calculation cannot 
be correct since for a single coupling theory the two loop /3-function is scheme independent. 

There are various well established ways of accounting for the effect of the evanescent oper- 



ators in the renormalization group functions. One approach developed in |2C] and used in the 
context of two dimensional four fermi theories is that of the projection or reduction formula. By 
examining the renormalized Lagrangian with the full set of evanescent operators incorporated 
one can relate any true renormalization group function to the naive ones which have evanes- 
cent effects present by adding terms in a well defined procedure which exactly accounts for the 
projected contribution from the evanescent operators, pO]. In essence the contribution of the 
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evanescent operator inserted in the appropriate Green's function but evaluated after operator 
insertion renormalization in strictly two dimensions, accounts for the overcounting in the nai've 
renormalization group functions. The alternative approach which we will use here and in |15| is 



to examine the renormalization group equation, which involves the true renormalization group 
functions, acting on the finite renormalized Green's function. As this is a purely two dimen- 
sional object the presence of the evanescent operators is effectively washed out and hence one 
can deduce the renormalization group functions by ensuring the renormalization group equation 
is valid to the order in perturbation theory one is interested in. However, for this approach to 
be successful one has to evaluate all the Feynman integrals at that loop order to the finite part 
inclusively. Ordinarily when renormalizing a theory at high loop order determining the finite 
part is a difficult aspect of the calculation. However, in using this strategy here it will turn 
out that it is possible to determine all the necessary parts of the four loop Feynman graphs 
to achieve this aim. Moreover, we have checked our final answers explicitly for the case G = 
SU{Nc) by using the projection formula of |2C] as an independent verification. 



Having reviewed the issue of evanescent operators in the d-dimensional extension of the 
NATM in relation to the construction of the renormalization group functions, we will now 
introduce their explicit forms for the calculation. First, we need to define the basis for the 
7-algebra in d-dimensions which has been introduced in |16, 18, 31 1. In d-dimensions the 
Clifford algebra 

{7^7"} = 27?^^^ (2.7) 

remains central to all 7-matrix manipulations except that the Lorentz indices now run over 
1 < /i < d where d is non-integer. This latter property means that, for example, the anti- 
symmetrized product of five or more 7-matrices is non-zero in d-dimensions whereas they would 
clearly vanish in four dimensions. Therefore, this set of this combination of 7-matrices is infinite 
dimensional in d-dimensions and can be used to span the space occupied by the d-dimensional 
7-matrices. For compactness we introduce the objects F^^j"^" which are anti-symmetric in the 
Lorentz indices and defined by 

pjn...M„ ^ ^[/ii _ ^ ^ ^Mn] (2.8) 

as the basis for the 7-matrices when we treat ( |2.1| ) in d-dimensions. Before considering a general 
group G it is instructive to recall the structure of evanescent operators for the case G = SU (Nc) 
as it will also illustrate other issues such as the restoration of multiplicative renormalizability. 
For SU{Nc) the group generators satisf 



rpa rpa ^ 

^IJ-'-KL — 2 



(2.9) 



This means that when performing loop calculations based on the interaction of (2.1) any string 
of group generators involving an even number of T'^'s can be written in terms of the SU{Nc) 
tensor basis of ( p.S| ) which is I I and T"" ^T"". Therefore, one can immediately write down 
the most general Lagrangian in d-dimensions which contains (|2.1| ) in two dimensions. We have. 



■ 7N 2 



2 -1 ^ 

^natm ^ i^^^ ^^^^ + I (^^p^i ^^^T^J^p^J) + - J] 5fco (V5*'r(,) W 

A:=0 

+ 2 E (2.10) 



^ k=0,k^l 



where all quantities are bare and the Roman letter k labels evanescent contributions in the same 
notation as |15]. As in [18, 15] we have introduced a coupling constant for each evanescent 



'Analogous decompositions exist for the other classical and exceptional groups. See, for example, 
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operator which thereby ensures ( |2.1[l| ) is now multiphcatively renormahzable in d-dimensions. 
Ordinarily one would use ( 2.10| ) to determine the true renormahzation group functions but from 



the point of view of calculabihty ( 2.10| ) is unpractical. To see this one need only consider the 



fact that at one loop each operator needs to be renormalized and given the nature of the one 
loop diagrams there will be mixing between all the operators so that proceeding to higher orders 
becomes a huge exercise. For instance, the vertex ['tpT'^^i^ijj)'^ generates the additional vertices 
{'ipIT(^^-^ip)'^ and {tpT'^T (^^•^ip)'^ at one loop, |17, |l5|. However, to achieve the aim of determining 
the two dimensional renormalization group functions this is not necessary since in the limit to 
two dimensions the new operators will be absent or equivalently one can set the new evanescent 
couplings to zero. In this case new operators will be generated under renormalization with 
a coupling which is the effective renormalization constant. As there will be a finite number 
of these at each new loop order it is much easier to calculate their effect and renormalization 
compared to an infinite set. Their full presence in the naive renormalization group equations are 
then accounted for by the formalism of [20, 17], discussed earlier. Having recalled the structure 



for the case G = SU{Nc) the general group situation is only complicated by the absence of a 
general identity (p.9|). In other words the complete set of operators analogous to those of ( 2.10| ) 



in SU{Nc) will still involve the r(„)-matrices but the group theory content has to be replaced 
by general functions of the group generators which are a basis for the group space. Whilst we 
do not know how to determine this to all orders in general we instead define the d-dimensional 
Lagrangian with evanescent four fermi operators as 

^natm ^ 01 ^ 9 f^U ^.^fj^^^J^^ + 9k a^Of (2.11) 

fc=o {a} m 



where 

O 



f = \ (i,^'T^,^T}^jr') (^^^■^r(fc)T^^V^-^) . (2.12) 

Here 7} j are functions of the group generators and the indices a and (3 represent sets of free colour 
group indices A. If one knew the full basis explicitly it would be possible to relate the SU (Nc) 
evanescent couplings gkap to those of the general case. As we will work with nullified evanescent 
couplings some insight into the possible form of the group basis can be gained by examining the 
renormalization of the four point function, |15]. This is in fact necessary here because we need to 



know which operators are generated as they have to be included in the subsequent order of the 
two point function diagrams and their contribution to the naive wave function renormalization 



accounted for. In [17| the new operators generated in the one loop four point function were 
calculated for general G and incorporated into the two loop calculation. However, the new 
operators generated at that order were not recorded. Whilst they were determined in at 
two and three loops for SU (Nc) these results were expressed in the {/ ® /, T""} basis which 
also cannot be uniquely generalized to arbitrary G. Therefore, we have recalculated the full 
two loop 4-point function renormalization for general G. This involved computing the graphs of 

figure 1 which generated the evanescent operator O32 = ^ {4'^^^{3)Tj'f^ip^'^^ at one loop where 
the operator is labelled according to the respective number of 7-matrices and group generators 
present. In figure 1 the vertex corresponds to the original vertex of ( |2.1| ) whilst for the full two 
loop calculation the additional graphs of figure 2 must be included where the symbol (X" at a 
vertex there represents the insertion of the operator ©32. We find that the full renormalized 
Lagrangian to two loops is 
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Figure 1: One and two loop corrections to the 4-point function. 



+ 



(2.13) 



where fl is the scale introduced to ensure that the coupling constant remains dimensionless in 
(i-dimensions[] and we note that the evanescent operators all involve the symmetric product of 
group generators except Z33. Here the fields and the coupling constant g are regarded as the 



N[032]]g 





SN[032] 



Figure 2: Single evanescent operator contributions to the 2 loop 4-point function. 



renormalized ones as opposed to those in {2A) and ( ^.lOD which were assumed to be bare. We 
find the wave function renormalization constant is 



1 



C2{R)T{R)Nfg-' 



C2{R) 



C2{G)T{R)Nf 
12e2 



+ 



[4r2(i?)iV/ + C2{G)T{R)Nf + 2Cl{G) - 6C2{G)C2{R) + 4C|(i?)] 

24^ 



and the four point operator renormalization constants to 0{g^) are 



1 + 



G2{G)g ^ ( Gl{G) _ [Gl{G)-SG2{G)T{R)Nf\ \ g^ 



27re 



^32 



>13 



53 



27re 

27r2e 
1 
4? 



+ 



4e2 

[7G2{G) - AG2{R)\ 



2e2 



+ 



32e J vr' 

7G2{G)-2G2{R)-2T{R)Nf] 



1 

16^ 



-^33 
7r2 



15 



15 



£2 + 2e 



51 



4e 

C7|(G) 



1 



192e^ 



^ + 0(5^) 
(2.14) 



+ 



384e 



(2.15) 



These values reduce to those quoted in [|T5[ for the case G = SU{Nc). Although to two loops we 
have managed to write the generated operators in terms of anti-symmetric or symmetric products 
of strings of group generators, it is not clear to us whether this should persist to all orders. For 
example, it may be the case that at a very large loop order the Lie algebra manipulations we used 
to obtain this Lagrangian do not allow us to reduce the generator products to pure T°-strings 

^We use jl as the renormalization scale throughout instead of the conventional fi since the latter will subse- 
quently denote d/2. 
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in that one or more structure functions may remain or three or more traces of generators. A 
systematic study which goes beyond the symmetric traces of group generators for general G of 
|27| would be needed. However, to the order we will be computing the wave function anomalous 
dimension, ( |2.13| ) is all that will be necessary if one considers the topologies of the diagrams 
which will arise and the operators which can be inserted at various vertices. 

We close this section by recalling the relationship of the NATM with other models. As is 
already well established the abelian Thirring model for = 1 is trivially equivalent to the 
single flavour Gross Neveu model due to the Fierz identity for the Dirac fermions we are using 
which is 

{i)-i>'i;f = - 2{iPipf . (2.16) 

Therefore, the renormalization group functions for both models in this limit must agree. How- 
ever, when one examines the MS renormalization group function in each case it transpires that 
beyond the leading order there is disagreement which arises because ( ^.16 ), which is established 
using the properties of the strictly two dimensional 7-algebra, is not valid in d-dimensions. 
Indeed one can use a d-dimensional Fierz lemma, ||3^ , 34|, to discover that an infinite set of 



additional operators emerge on the right side of (2.16) in this case. Therefore, to restore the 
clear equivalence between the renormalization group functions in both models one needs to make 
a finite scheme change in the NATM results. Introducing a finite renormalization to ensure a 
symmetry principle is preserved in the quantum theory is not novel to the NATM. For example, 
the renormalization of the flavour singlet axial vector current in QCD in MS ceases to preserve 
the chiral anomaly and one ensures full consistency by a finite renormalization, Whilst 
it may appear here that the problem is not related to the 7^ issue of the QCD axial anomaly, it 
is in fact implicit in the computation since the d-dimensional object T'^^) would be the projector 



of the two dimensional 7^. Therefore, as in |l|] it will be the case that having established an 
MS result for the four loop anomalous dimension we will need to determine the constraints on a 
class of possible scheme changes to preserve the Gross Neveu equivalence. In [jl^ another equiv- 
alence was mentioned between the Nf = 3 Gross Neveu model and the = 1 5'C/(4) NATM. 
Whilst the implications of this were discussed for the three loop NATM /3-function, we will 
not impose it or discuss it here primarily because it is not clear if a full rigorous proof of the 
relationship has been established. If one does exist it will be a trivial exercise to determine its 
consequences for the anomalous dimension. For later use and to illustrate these remarks we now 
quote the relevant renormalization group functions in MS in both models to the orders they are 
known. For the Gross Neveu model the anomalous dimension in our notation is |22, p3|, 24], 



7(5) 



and the /3-function is, pll. 



mf - l)g2 {2Nf - l){Nf - l)g^ 



+ 



{2Nf - 1){4N^ 



167r3 
14A^/ + 7)g^ 



1287r4 



+ o(/ 



Pig) 



{Nf - l)g^ ^ {Nf-l)g' ^ {2Nj - 7){Nf - l)g^ ^ 



vr 27r2 IGvr^ 

where we note that the Gross Neveu Lagrangian is, pl| , 

• \ 2 



L 



GN 



(2.17) 



(2.18) 



(2.19) 



For the NATM the anomalous dimension is, [14, 15, 18, 19 
C2iR)T{R)Nfg^ 



7(5) 



27r2 
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2C2{G)C2{R) - CUG) - 2C2{G)T{R)Nf - 8T\R)N/ 



167r3 



+ 0{g' 



(2.20) 



and the /3-function is, fl], |T|, [T|] , 

C2{G)g^ ^ T{R)NfG2{G)g^ 

"5 



27r 
+ G2{G) 



27r2 



T\R)Nf + ||c|(i?) 



67 



31 



^2C2(i?)C2(G) + f^CliG) 



(2.21) 



We also note that several all orders results are available in the Thirring model, though it is not 
clear if they are in an MS scheme. In the abelian case, |3£], 



7(5) 



Nfg^ 



2vr2 (1 - ^ 



Pig) 







(2.22) 



For the NATM, Kutasov has written down an expression for (3{g) which has a similar form to 
(|2.22| ), [^ ]. In the current notation of this paper it is 



Pig) 



C2iG)g^ 



2tx (] 



TiR)Nfg 



2lT 



(2.23) 



This expression was established through a connection the two dimensional theory has with string 
and conformal field theory but was in fact only valid when Nj is large. More recently, however, 
by using current algebra Ward identities this result has been argued to be exact at all orders 
in a particular scheme, [^l[. This is an important point since it only contains the elementary 
Casimirs of G and not the higher order ones which we are concerned with here. Whilst ( 2.23| ) 
would put additional constraints on potential scheme changes needed for the /3-function we will 
not need to consider these here as they will not play a role in the analysis of the anomalous 
dimension. 



3 Four loop calculations. 

We are now in a position to carry out the four loop renormalization of the wave function in 
(2.1) using the d-dimensional evanescent Lagrangian of ( 2.13| ). First, we redo the three loop 



calculation of |T|, 0, |T8|, |l|l but for the case of a general Lie group. Unlike ||l5| we will use the 
massless propagator 



for the fermion field. In |15] the massive NATM was renormalized to three loops using a fully 
massive propagator and a propagator which was massless but which included a mass term in 
the denominator to act as an infrared regulator. This was necessary since that paper considered 
the renormalization of the full model including the four point function and not just the two 
point function as we are here. For the former Green's function such an infrared regulator is 
necessary to avoid potential spurious infrared infinities which could arise when using completely 
massless propagators with two of the external momenta nullified. In the presence of a mass with 
all external momenta nullified means that the poles in e which emerged represented ultraviolet 
infinities only. For the two point function in the fermionic theory we are considering such 
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spurious infrared infinities ought not to arise at the order we are working to since, for instance, 
neither external momentum is nulhfied. Therefore, in all our calculations we use (^]^) without 
an infrared mass regularization and this allows us in fact to calculate all the necessary integrals 
exactly to four loops or at least to the finite part in e. This latter property is important since 
we will be considering scheme changes later and it is necessary to have the renormalized Green's 
function in order to extract the renormalization group function in various schemes. Further, 
by choosing a massless propagator the Feynman diagrams with tadpoles| are absent since by 
Lorentz invariance the tadpole subintegral is zero. Therefore to three loops the only two basic 
topologies which are relevant for the wave function renormalization are those given in figure 3. 
They can be evaluated exactly as a function of d = 2/x and we find that they are respectively 




Figure 3: Two and three loop corrections to the 2-point function. 



_ iC2{RmC2iR) - C2iG))ifi - 1) + 2TiR)Nf]T{3 - 2^i)T\^,)T{^, + l)g^ 

(/x-l)r(3/.-l) 

and 

16 

8{3i2^ -9fi + A)C2{R)C2{G) - — (V - + 5)C|(i?) - 16{fi - l)C2{R)T{R)Nf 
r(4 - 3^)r2(2 - ^)r(3/i - 2)THi,)C2{R)g' 

{^i - i)r(3 - 2/i)r(4^ - 2)r2(2^ - 1) ^ ' 

where we have omitted here and in later values of Feynman graphs, the overall power of the 
momentum and renormalization scale fx, which can be restored on dimensional grounds, |5, where 
p is the external momentum, and factors of (47r)°'/^. Also the coupling constant g which will 
appear in the values is bare. However, for the full three loop anomalous dimension to be deduced 
correctly we need to take account of the contribution from the generated evanescent operator 
032 which is represented by the Feynman diagrams of figure 4. In these graphs the 7-algebra 
now involves the F (3) -matrix sandwiched between ordinary 7-matrices. To handle this in the 





Figure 4: Single evanescent operator contributions to the 3 loop 2-point function. 

computation we decompose r(3) into its six terms, multiply the 7-string by ^ and take the 
trace as well as dividing by the normalization of 2p^. To cope with the algebra generated and 
the mundane massless integration in this and other graphs we have made use of the symbolic 

"''By tadpole in a theory with a four point interaction we mean diagrams with zero momentum insertions. 
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manipulation programme FORM, [25|. Moreover, to ensure that the integration routines are 
properly constructed we have repeated the full four loop anomalous dimension calculation of 
p2| in the SU{N) Gross Neveu model and reproduced the result of |22]. Having verified the 
programmes reproduce the correct results for this model it is then elementary to re-run them 
where the input Lagrangian is replaced by the NATM. Whilst evanescent operators will not be 
present in the four loop Gross Neveu calculation the integration routines were still valid for the 
NATM since we evaluated systematically all the basic scalar integrals which could arise in any 
of the topologies. Although the decomposition of r(3) for the graph of figure 3 was relatively 
efficient, for higher order F^^) 's or for diagrams with more than one evanescent operator insertion, 
it was much more appropriate to make use of the general properties of the F(„)'s in these cases. 
For example, the lemma, [17, 18, 31], 



(n) (m) Mi---Ain 



/(n,m)F^^)-'^™ 



where 



Qn 



f{n,m) = (_l)"-(_l)n(n-l)/2 _ + 



(3.4) 



(3.5) 



u=0 



led to the other graphs being determined relatively quickly. Hence, we record the value of the 
graphs of figure 4 are 



iC2{R)[3C2{G) - 4C2iR)][C2iG) - 2C2(i?)](2^ - 1)(/^ - 3)/ir(3 - 2f,)T^fi)g^ 

F(3^ - 1) 



(3.6) 



We now turn to the details of the four loop calculation. The basic topologies are given in 
figure 5 and fall into two classes. The upper two Feynman diagrams correspond to elementary 







Figure 5: Four loop corrections to the 2-point function. 



massless chain integrals and can therefore be computed straightforwardly. We find their d- 
dimensional values are respectively 



3iC^{R) 



2C2{G)C2{R) - 2Cl{R) 



AC2{R)T{R)Nf 



1 



Cl{G) 



+ 



2G2{G)T{R)Nf 2T'^{R)Nf 



^^F(5-4/i)F5(^)g4 
(3;u - 2)F(5/i - 3) 



(3.7) 
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and 



12(32^^ - 213^^ + 581;u^ - 514^ + 144)C2 (G)C|(i?) 
4(56/i^ - 301;U^ + 701^^ - 586/i + 160)C|(.R) 

- 96{3n - 2){fi - l)Cl{R)T{R)Nf + 32{3fi - 2)C2{G)C2{R)T^iR)N^ 
(232/ - 1903^^ + 5791^^ - 5294/i + 1504)C|(G)C|(i?) 

+ 96i3fi-2){fi-l)C2iG)Cl{R)T{R)Nf 



16{3n - 2)C2{R)T^{R)N^ 



+ 



(296/ - 2903/x^ + 9479^^ - 8830/x + 2528)C|(G)C2(i?) 

6(/x-l) 

(432/ - 2004/ + 3863/ - 3693/ + 1734/i - 320)Cl{G)C2{R)T{R)Nf 

6(/x- 1)3 

4(8/ - 125/ + 461;u^ - 442// + 128) d'^'^^df'"^ 
- l)iVfund 

2(12/ + 119/^3 - 237/^2 + 150/i - 32)df'"^df'"^Nf 
(/X - l)3iVfund 

r(5 - 4fi)T{Afi - 3)r3(2 - /r7(// 



64(3// - 2)Cl{R)T'^{R)Nf 



{2fi - l)2r(5/u - 3)r(4 - 3/r3(2/u - 1) 



(3.8) 



Although these are elementary to derive we draw attention to the group theory structure of 
each result. The former since it has a propagator correction involves the usual products of 
Casimirs which arise at lower order. By contrast the embedded chain diagram involves the two 
new Casimirs constructed from the tensors dp and dA- It is instructive to consider how these 
arise. Clearly the term df^'^df'"^ arises from two strings of group generators where one of the 
T'^ from T°'®T"' tensor structure is located in one string with the other in the second string. On 
the other hand the second new combination arises from manipulating particular strings of eight 
generators. They are ah related to Tr (^arpbrpcrpdrpaj^hrpcrpd^^ Using the Lie algebra, (|]2|), it is 
elementary to deduce 



rp^^ / rparj-ibrj-iCrpdrparpbrpCrpd 



.r=0 



Y[(G2iR)-"-C2{G)) - cl{G)C2{R) 



^fund 



_j_ japq jbqr jcrs jdsprj^^ f rparpbrpcrp, 



(3.9) 



Examining the last term reveals that the combination of structure constants is equivalent to 
Tr (a^A^^M^) where A" is the adjoint representation of the generators and hence this string 

can be related to the d'p'^'^d'^'^'^ term. One can also consider the group theory of this diagram 
from the point of view of the connection with QCD. In (|2.1| ) it is possible to write the four 
point interaction in terms of a spin-1 auxiliary field which in d-dimensions is related to the 
gluon of QCD. In such a reformulation the second diagram of figure 5 would correspond to 
the first graph of figure 6 where the spring line is the auxiliary field. Clearly one can identify 
the d'p'^'^d'p'^'^ structure that is contained in the topology. However, it is not evident how the 
^ahcd^abcd g^ructure would correspond. For this case one would have to have a diagram akin 

to the second graph of figure 6 where the gluonic loop would lead to Tr [A'^A^A'^A'^^ . As this 
diagram cannot be reduced to any of figure 5 by restoring the original interaction it might appear 



13 



Figure 6: Four loop self-energy corrections in QCD containing d'p^'^ d'p'^'^ and d'p'^'^ d°^^'^ respec- 
tively. 

that there is a potential inconsistency in our evaluation. However, as will become apparent later 
this issue will be satisfactorily resolved. For the remaining two diagrams of figure 5 we expect 
^abcd^abcd ^abcd^a^cd iqj-t^^ present after considering the routing of the group generator 

strings in the graphs. 

The main point concerning these two diagrams is that their evaluation rests in the fact 
that unlike the previous graphs discussed so far they do not reduce to simple chain integrals. 
Instead upon taking the spinor trace and evaluating two elementary loop integrations one is 
left with a set of two loop integrals which need to be evaluated. These have the general form 
of figure 7 where we use coordinate space representation where the internal momenta k and I 



k 




Figure 7: General two loop self energy topology whose value is denoted by (ai, 02, as, 04, as). 

are integrated over. The ai beside a line denotes the power to which that scalar propagator 
is raised to. The graphs either involve tensor numerators or are purely scalar. For the former 
case given the large symmetry of the diagram one can always rewrite such numerator scalar 
products in terms of denominator factors which then reduce the exponent of that line by unity. 
After this procedure one is left with diagrams which are either chains or purely scalar two loop 
graphs. For each of the graphs of figure 5 this remaining topology has the same pattern of 
exponents. For the third graph we find {ai} = 1} and for the final graph {ai} = 

{1, 1, 1, l,a} where a and (3 here are functions of d and therefore non-integer when d = 2 — e. 
By integration by parts rules and recurrence relations based on the uniqueness rule of ^] 
one can reduce this set of graphs either to chain diagrams or to two respective d-dimensional 
basis scalar integrals given by (1, 1 — /i, 1, 1 — /x, 1) and (1, 1, 1,1,1— fi) where we have again set 
d = 2/i for shorthand. These then are the outstanding diagrams which need to be determined. 
Given the form of the integrals it turns out that they cannot be evaluated in closed form for 
arbitrary values of d. However, their e expansion can be deduced to an order beyond that which 
we need for this calculation. This is achieved by using the Gram determinant method of [^] 
where the value of a d-dimensional two loop self energy diagram of the form of that of figure 7 
can be exactly related to a ((i+2)-dimensional diagram of the same form plus a set of elementary 
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chain integrals. In [^] the Gram method was applied to a particular choice of the {ai}. In the 
appendix we briefly recall the arguments and record the general expression for the relation. In, 
for example, the application of the rule of (|A.4|) to the set of the form (1, 1, 1, l,a) the ninth 



and tenth terms of (A^) correspond to the same original topology. The actual combination of 



exponents can be restored by the recurrence relations of [44|. Having related these two basic 
integrals to ones in (4 — e) -dimensions, the e expansion of the higher dimensional integral is 



determined from the results of [45, 46 1. Whilst the e-expansion is known for the general integral 
to O(e^), [45 1, only the first term is relevant to our two dimensional calculation as the factors of 
— 1) ensure the leading term appears in the finite part. It might be thought that the group 
theory and symmetry arguments used in [^, |4^] to deduce the e-expansion of the problematic 
two dimensional integrals could be used directly on them. However, these arguments seem to 
rely in part on the fact that in four dimensions the diagram is finite with respect to e and 
the leading term of the expansion is independent of {ai} for {ai} of 0(e) from unity. In two 
dimensions the diagram is singular in e and as the residue of the pole depends on the parameters 
it does not appear possible for the group theory arguments of ^ to proceed. Given these 
remarks we record the value of the third graph of figure 5 is 



(8(49/ 



336/ + 794/ - 832^ + 283)C2(G)C|(i?) 



- 8(27/ - 174/ + 392^2 _ 394^ _^ 131)G|(i?) 

- 8(^- 1)(18/ -59;U + 32)C|(i?)T(i?)iV/ - 4(6/ - 15// + 8)G|(i?)r2(i?)A^/ 

- 4(59/ - 430/ + 1072/ - lUSfi + 411)C|(G')C'|(i?) 
+ 4(/i - 1)(42/ - Wlfx + 88)C2{G)CUR)T(R)Nf 

+ -(36/ - 281/ + 747/ - 867/x + 311)G|(G)G2(i?) 
3 

(576/ - 3604/ + 6853/ - 5159// + 1346)G|(G)G2(i?)T(i?)iV/ 

12(/x-l) 

+ 2(6/ - 29fi + lQ)C2{G)C2{R)T'^{R)Nf 
8(3/ - 38/ + 144/ - 210/i + 82,)df'"^df'"^ 
^fund 

2(4/ - 37/ + 23// - 2)df^'^df^'^Nf\ r^{3 - 2//)r6(/x) 



(/. - l)iVfund J (P- 1)2(2/^ - l)(5/x - 4)r2(3/x - 2) 

+ (384(77/ - 605/ + 1743/ - 2395/ + 1498// - 346)(// - l)G2(G)G|(i?) 

- 384(41/ - 304/ + 842/ - 1124//^ + 691/i - 158)(// - l)G|(i?) 

- 768(12/ - 47/ + 49/i - 16)(// - lfCl{R)T{R)Nf 

- 192(97/ - 804/ + 2414/ - 3424/i2 + 2183// - 510)(/i - 1)C^{G)C^{R) 
+ 1536(7/ - 32/ + 34/i - 11)(// - lfC2{G)Cl{R)T{R)N} 

- 1536(// - lfGl{R)T'^{R)Nf 

+ 32(126/ - 1103/ + 3467/ - 5093/ + 3315/i - 784)(/i - l)G|(G)G2(i?) 

- 8(384/ - 2713/ + 6378/ - 6806//^ + 3431// - Q7Q)Gl{G)G2{R)T{R)Nf 
+ 192(4/ - 23/ + 25// - 8)(// - l)G2{G)G2{R)T'^{R)Nf 

Mbcd^abcd 

- 384(9/ - 100/ + 382/ - U<^^? + 447// - 110) (// - 1) \^ ^ 

^fund 

.4 , 2 ,df'"^df'"^Nf\ r2(2-/r4(//) 

- 192(/-18/ + 14/ + /i-2)--^ ^ n V VP; 



Niund ; 6(2//-l)(5//-4)r2(2//-l) 
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X (1,2-^,1, 2-/x,i; 



M+1 



+ 



/ 4(19^4 - 126/i3 + 294^2 _ 394^ + 105)C|(i?) 



+ 16(3/^2 - 10^ + 6)Ci{R)T{R)Nf 



+ 



+ 



4(35/i^ - 247^^ + 603/i^ - 647^ + 228)C2(G)C|(-R) 
2(43/i^ - 322^^ + 826/i^ - 924^ + 333)C|(G)C|(i?) 

4(14^2 _ 55^ ^ Z2>)C2{G)Cl{R)T{R)Nf + 4(2/^ - 3)C^{R)T^{R)N^ 

(54//^ - 431;u^ + 1171/i^ - 1373// + 507)g|(G)C2(-R) 
3(/x - 1) 

(192/ - 1225/ + 2387/ - 1854^ + 50A)C^{G)C2iR)T{R)Nf 

12(/i-l)2 

4(/ - 5^ + 3)C2{G)C2{R)T\R)Nf 



+ 



+ 



4(3/ - 34/ + 122/x2 _ 172^ + 69)df'"^df'''^ 
in - l)iVfund 

ifX - l)2iVfund j (M-l)(2^-l)r(5^-3) 



(3.10) 



The term (1, 2 — 1, 2 — ^, 1)|^_,_]^ corresponds to the higher dimensional integral derived in the 
Gram determinant method where the restriction denotes a 2(// + l)-dimensional measure. Space 
prevents us from recording the exact value of the remaining graph since it is of a comparable 
size to (|3lO| ) though we note it has a similar form. 

To complete the full four loop calculation the contributions from the evanescent operators 
need to be included. One can divide these into the diagrams which involve the operator generated 
at one loop in the four point function and those arising from the two loop ones. For the latter, 
the relevant topologies are those of figure 2. We record the respective values of the graphs and 
include the appropriate renormalization constant to aid identification. They are, in addition to 
which involves Z32, 



iZ 
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9C2{G)CUR) - 4CliR) 



on , 

-CI{G)CUR) + -Cl{G)C2{R) 



2^a,bcd ^abcd 
^fund 



5^ 

^r(3-2/r3(// 



(/i - l)iVfund 



r(3^ - 1) 



4iZ3iC2iR)[2C2iR) - G2iG)] 



/2^-l)(/i-3)r(3-2/r3(// 

r(3^ - 1) 



(3.11) 



(3.12) 



33 



';C2{G)CliR) 



lciiG)Gi{R) + I 



+ 



^(2/x - 1)(^ 



Gl{G)G2iR) 

3)r(3 - 2^)r3(^)/ 



r(3^ - 1) 



8iZMR)[2G2{R) - C2(G)]/i(2/i - 1)(m - 2)(^ - 5) 



r(4-2/r3(// 

r(3/. - 1) 



(3.13) 
(3.14) 



16 



and 



53 



18C2(G)C|(i?) - 8C|(i?) 

M(2//-l)(M-5)r(5-2^)r-^(^)5- 



^Cl{G)Cl{R) + ^C|(G)C2(i?) 



-^fund 



r(3/x-i) 



(3.15) 



These explicit values illustrate that even though they are two loop topologies some do contain the 
new Casimir structures that appear at four loops. This is because unlike the original interaction 
these operator insertions have more than one pair of group generators. Hence the strings of 




Figure 8: Double evanescent operator correction to the 4 loop 2-point function. 

group generators can contain one set of eight T"'s. The remaining set of contributions come 
from the insertion of O32 in the topologies illustrated in figures 8 and 9. For the three loop 






Figure 9: Single evanescent operator correction to the 4 loop 2-point function, 
topologies their values are 

"32(4/; - r){^i - SfC^iR) 8(34/,2 - I8i/i + 205)(/; - 3)C2(G)C|(i?) 



(m-1) 

+ 32{ii - 3)Cl{R)T{R)Nf + 



- 40{^i-3)C2{G)Ci{R)T{R)Nf + 



2(73/^2 - 462/x + 545)(// - 3)C|(G)C2(i?) 
2{6n^ - 24/x + 19)C^{G)C2{R)T{R)Nf 



4(49/^2 - 287/x + 332) (/X - 3)C|(G)C|(i^) 



(M - l)^fund 



+ 



+ 



^Qdf'^df'^Nf 
(m - l)-^fund 



(2/i - l)r(4 - ?,ii)T\2 - /i)r(3^ - 2)V\ii)g^ 



r(3 - 2iu)r(4/[i - 2)r2(2/i - 1) 

Whilst the diagram with a double operator insertion has the value 

20(2/^2 - 15^ + 19)(Af - 3)C2(G)C|(i?) - 8(2^^ _ ^ -^^-^^^ _ ^)C*{R) 
33(2//2 _ 15^ + 19) (^ _ 3)C|(G)C|(i?) C'i{G)C2{R)T{R)Nf 



(3.16) 
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+ 



55(2/i2 _ 15^ + 19) _ 3)C'i{G)C2{R) 8(2/i2 - 15/i + 19) 



3 ) dP^'^'^ ^ahcd 



fund 



12 

/i(2/i 



fund 



i)r(3 - 2/i)r3(^)52 



r(3M - 1) 



(3.17) 



This diagram contributes when one considers the powers of the couphng constant which are 
present. This completes the evaluation of all the diagrams relevant for the renormalization of 
the two point function at four loops. 



4 Four loop anomalous dimension. 



Having determined all the contributions to the two point function including those from evanes- 
cent operator insertions it is a relatively straightforward task to determine the four loop correc- 
tion to the wave function renormalization. We find 



.MS 



1 



4^e ^'^^> [ 

+ [AT\R)N^ + C2{G)T{R)Nf + 2C|(G) - 6G2{R)G2{G) + 4C|(i?)] 



24e 



TiR)NfG2iR)GliG) 



32e3 



+ 



192- 



^abcd^abcd 
^fund 

+ 352GUG)G'i{R) + 3Gl{G)G2{R)T{R)Nf 
+ 32G2{G)G2{R)T'^{R)Nf + 192C7|(i?) 



48 



Nfdf^'^df'"^ 
^fund 

USG2{G)Gl{R) 



96GliG)G2{R) 



8Gi{R)T^{R)N^ 



256e^ 



(^mbcd Jabcd jabcd jabcd at 

772Cl{G)G2{R) - 1536%^ - 624^^-^^^ 
^fund ^fund 

- 2808C| (G')C|(i?) + 2mGl{G)C2{R)T{R)Nf + 3552C2(G)C|(i?) 

- SQ^G2{G)Gl{R)T{R)Nf - 1536C|(i?) + 57QGl{R)T{R)Nf 



+ 192C2(i?)r3(-R)iV/ 



1 



1536e 



+ 0(/) • 



(4.1) 



Ordinarily one does not quote this value but merely gives the corresponding anomalous dimen- 
sion. We have done so here to illustrate an important feature of renormalizing a model with 
evanescent operators. For ordinary (single coupling) theories where this is not an issue the poles 
of the MS renormalization constants are related in a particular way. For example, the residues 



of the non-simple poles are determined by the simple pole residues at lower orders, |47|. How- 



ever, examining (4.1) one can see that this is not the case since at four loops the new Casimir 
contributions occur in the residue of the double pole in e. Usually this would indicate an error in 
the renormalization since these new Casimirs cannot arise before four loops and so they should 



only appear in the simple pole in e. In this instance, however, the choice (4T) correctly renders 
the two point function finite but the anomalous dimension derived from it corresponds to the 
nai've renormalization group function and not the true one in relation to the discussion of section 
2. To obtain the true anomalous dimension there are two possible procedures to follow. One 
is to extend the projection technique developed in pO, |l^ and applied to the SU{Nc) NATM 



/3- function in |15] to the present case. This would involve calculating the corrections to the pro- 
jection functions by inserting each evanescent operator into a two point function and calculating 
the finite part after renormalization. These finite parts determine the projection function for 
each operator and together with the naive /3-functions of each operator, allows one to deduce the 
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true renormalization group function. The alternative approach which produces the equivalent 
result is to ensure that the finite renormahzed two point function in two dimensions satisfies the 
renormalization group equation 



d d 71 







(4.2) 



where T^"'\p, fi,g) is the renormahzed n-point Green's function. To ensure this condition can 
be used one must be careful in retaining the finite parts of all the integrals in the construction 
of the Green's function before renormalization. This is the reason why we have been careful to 
compute the integrals exactly in most cases and to sufficient powers in e in the other cases and 
means we will follow the latter course here. After renormalization the two point function is 
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^fund 
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1U7GUG)G2{R) + 4452C7|(G)C|(i?) + n52C2{R)T^{R)Nf^ 

6096C72(G)C|(i?) 



- 81^Cl{G)C2{R)T{R)Nf 
+ U08G2{G)Cl{R)T{R)Nf 
+ 2856C|(i?) - 2688Ci{R)T{R)Nf 



62AC2{G)G2{R)T^{R)N? 



+ (576C(3) 



3312)-^—^ ^ 

^fund 
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30727r4 



lQbQCl{R)T^{R)Nf 
4 + 0(/) 



(4.3) 



where /i^ = 47re and 7 is the Euler-Mascheroni constant. Hence, using the three loop MS 
/3-function of p^] we find that for a general colour group the MS result is 



7^^(^7) 



C2{R)T{R)^^ 



+ 



2C2{G)C2{R) - GliG) - 2G2{G)T{R)Nf - 8T'^{R)Nf 



624 



^fund 



+ 57C|(G)C2(i?) - 198C|(G)C|(i?) 



C2(i^)g' 
167r3 



8'iGl{G)C2{R)T{R)Nf + 168C2(G) C|(i?) + UAC2{G)Gl{R)T{R)Nf 
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- 192C2(i?)r3(i?)iV/] ^ + 0(/) . (4.4) 

As a check on this value we have in fact applied the projection formalism to the case when 
G = SU{Nc) and verified that both are in agreement. One interesting feature of the expression 
([4.4| ) is that the new Casimir d'^'^'^d'^'^'^ has cancelled in the final expression leaving only a term 
involving d'p'^'^d'p'^'^. This cancellation appears to be consistent with our earlier discussion on 
the nature of the graphs in the two point function when regarded from a QCD point of view and 
also reflects the contribution from the relevant evanescent operators. Though we do not regard 
this observation as a hard check on the final result. More importantly from the point of view 
of our original motivation a new Casimir does appear in the (MS) anomalous dimension of the 
fermion. 

Whilst we have produced the MS anomalous dimension at four loops there remains one final 
task to perform which has been discussed in related work. This rests in the nature of the dimen- 
sional regularization. In two dimensions one has various equivalences between certain models 
such as the abelian Thirring model being the same as the Gross Neveu model for Nj = 1. This is 
established through an elementary two dimensional Fierz transformation. In d-dimensions, how- 
ever, such relations cannot be preserved. In the first instance, the Fierz transformation becomes 
infinite dimensional as one has to decompose I I into the full basis of tensor products of T^^^, 



1 33, 34|. Secondly, the relation that one can establish depends on the evanescent operators which 
are generated in the four point interaction but not in a way which the direct relationship can be 
determined. Therefore the situation is such that in renormalizing the NATM, the resulting MS 
renormalization group functions do not preserve the specific equivalences. In other words taking 



the abelian limit of ( |4.4] ) does not recover the four loop result of [22| for the Nf = 1 Gross Neveu 
model. To ensure that this property is preserved in the renormalization having been broken by 
the regularization, we need to make a finite scheme change. Such changes were considered in all 
previous work in this area, [16, 17, 18, and we proceed along similar lines to |15| here. 

As the choice of possible scheme changes is infinite we restrict our study to the class intro- 
duced in |]l5|. There the scheme was changed from MS to one where the finite part was also 
absorbed into the renormalization constants. In order to examine the equivalences the finite 
part is parametrised and constraints placed on the parameters by ensuring agreement of the 
renormalization group functions in the various limits. Whilst this may introduce a large degree 
of redundancy it allows for choices in future applications of the results. Here as we are dealing 
with a general group G we use the usual Casimirs as the basis for the parametrization of the 
finite part. In particular we choose the finite part of to be 

C2{R) [W21C2{R) + W22C2{G) + W23T{R)Nf] ^ 

+ C2{R) [w3lGl{R) + W32Gl{G) + W33T\R)Nf + W3iC2{G)C2{R) 

+ W35G2{R)T{R)Nf + W3QC2{G)T{R)Nf 



4 + 0{g') (4.5) 



where the {w2n} and {w^n} are constants which can be easily related to the SU{Ni.) choice of 
It is worth noting that choosing these values one has to redo the renormalization systematically 
at each loop order as the e-pole structure of the new becomes dependent on the parameters. 
To determine the new renormalization group function in this general scheme one can either follow 



the general procedure for scheme changes given in, for example, |28| which involves calculating 
the function which relates the coupling constants in both schemes or apply the method used to 
deduce ( [4. 41) from the renormalization group equation itself. In fact we used both methods to 



ensure that our manipulations are consistent. Hence, we find that for the general scheme the 
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anomalous dimension is 

jS'-(-g) = - C2{R)T{R)Nf^ 



C2{R) [{hi + i) T'iR)Nf + h3C2iR)T{R)Nf + {w23 + bu + i) C2iG)T{R)Nf 

3 



+ {W21 - i) C2{G)C2{R) + {W22 + 1^) Cl{G) 



yj-3 



\ (hl^ + 36ii + 2621 + 1) C2{R)T^{R)Nj 
+ {b23 - 2W23 + i^is + 611613) Cl{R)T\R)Nf 

+ (611612 + |6ii + |6i2 + 622 + W23 + IW33) C2{G)C2{R)T\R)Nf 

+ (626 - 2«;2i + ^6^3) Gl{R)T{R)Nf 

- {2w22 - iu;35 + 1-W21- 625 - §613 - 612613 + §611) C2{G)Gl{R)T{R)Nf 

+ (l«'36 + §^ + W22 + 624 + i6?2 + f 612 + i|6ii) Gl{G)G2[R)T{R)Nf 

+ (1^31 - ^ - §613) C2{G)Cl{R) + (i«;34 + i + i|6i3 - §612) C|(G)C|(i?) 

' '^W32 + il&12) Cl{G)C2{R) - l^2r.A 



4 + 0{f) (4.6) 



+ 



where g denotes the coupling constant of the new scheme and we have taken the finite part of 
the renormalization constant Z'^Zg to be 

[bnT{R)Nf + 6i2C2(G) + 6i3C2(i?)] - 

vr 

b2iT'^{R)Nf + b22T{R)NfG2{G) + b23G2{R)T{R)Nf 

+ 624C|(G) + 625C2(i?)C2(G) + 626C|(i?)l 4 + 0(5') • (4.7) 

Whilst we have discussed the constraints on the three loop parameters in we will redo 
that calculation here. This is partly to do with the fact that there were restrictions from two 
equivalences and for reasons already discussed we will only consider the one relating the abelian 
Thirring model and the Gross Neveu model. For this we recall that the abelian limit of the 
Casimirs is 

C2{R) ^ 1 , C2(G) ^ , T{R) 1 , df^'^df"'^ 1 . (4.8) 
Therefore, we find the constraints 

611 + 613 = - ^ 
621 + 623 + 626 - 2w2i - 2w23 = ^ (4.9) 

o 

Whilst these do not restrict the values of the parameters substantially any (numerical) choice 
must satisfy these equations for the two dimensional symmetry broken in d-dimensions to be 



established. The first equation of (|4.9| ) agrees with that of [15| when the parameters are converted 
to the case of G = SU{Nc). 



5 Discussion. 

One of our initial motivations was to examine the appearance or otherwise of new Lie group 
Casimirs in the NATM fermion anomalous dimension which are known to arise at four loops 
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in QCD. Whilst one of the two potential new Casimirs is absent in (|4.4| ) the one involving the 
fundamental representation remains. Therefore if the connection between the WZWN model 
and the NATM critical exponents is valid then the problem of how such a term cancels in the 
final critical exponent when approached from the renormalization group equation point of view 
is still open. One possibility is that of choosing a renormalization scheme in such a way that 
this new term is absent from (4^). Although this would appear to resolve the issue it has 
implications for the other renormalization group functions as this scheme choice impinges on 
how one renormalizes their associated Green's function. In particular this will impact upon the 
/3-function and it is possible they will be transformed into it. Even if this were not the case, the 
computation of P{g) itself at four loops in the NATM will produce these new Casimirs as well, 
which can be seen from examining the group theory of the Feynman diagrams of the fermion four 
point function. This would appear to contradict the recent exact /3-function of |41]. However, 
in that case it seems that the use of conformal symmetry has somehow excluded the higher 
order Casimirs. A similar feature occurs in four dimensional gauge theories at four loops. If 
one examines the four loop MS /3-function of QCD, |12], for an arbitrary colour group the terms 
involving dp and cIa will vanish in certain cases. For instance, when the quark is in the adjoint 
representation and Nf = ^ then the resulting /^-function contains no higher order Casimirs. This 
is, of course, the result of the theory possessing a new symmetry which is A/" = 1 supersymmetry. 
Thus, in the NATM it may be that a similar mechanism such as the two dimensional conformal 
symmetry used in the construction of [41| is responsible and powerful enough to exclude the 
^c^cd^abcd ^g^j^ ^i^g NATM /3-function to all orders in a particular renormalization scheme. 
Though it is not clear in this approach what form the all orders fermion anomalous dimension 
would take. Alternatively from the critical exponent point of view since P{g) will contain 
information on the non-trivial fixed point, Qc, at which the renormalization group functions 
are evaluated at to obtain the critical exponents, to examine the mechanics of the Casimir 
cancellation further would require the explicit form of the four loop NATM /3-function in order 
to carry out a test of this point of view at this order of approximation. The complexity of such 
a calculation is on a footing equal to the renormalization of the two point function at five loops. 
Moreover, the four loop /3-function of the usual simple Gross Neveu model, which would serve 
as a preliminary to a similar calculation in the NATM, has yet to be performed. Therefore, to 
understand this further would require a substantial amount of new calculations beyond those 
performed here. 
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A Gram determinant. 

In this appendix we briefly recall the application of the Gram determinant method of to 
relate a two loop self energy integral in ((i-|-2)-dimensions to a set of d-dimensional integrals. The 
method relies on several properties. First, the d-dimensional measure of a (massless) Feynman 
integral can be related through 

d^+^x = ^d'^x . (A.l) 
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Next the Gram determinant of three Lorentz vectors, x, y and z, is defined by 



Gr(x,y,z) = det 



/ xy xz 
yx y"^ yz 
y zx zy z^ 



(A.2) 



and through considering the Lorentz invariance of each integration, one can show that, |22|, 



y" dS^x J d'^y J d'^z GT{x,y, z)F{bix^ + b2y^ + bsz^ 



d{d-l){d-2) 



/ d'^+^y / + ftay^ + hz^) 



(A.3) 



where F{x) is a general function. In our case it represents the Feynman parametized two loop 
self energy graph of figure 7 when converted into a vacuum diagram. For this topology one 
can always write the integral as a linear combination of the three Lorentz invariants x^, and 
where y and z represent the internal k and / momentum integrations and x corresponds to 
an integration over the endpoint to produce the three loop vacuum bubble. Since the Gram 
determinant is invariant under the change of variables used to produce the form of the argument 
of F{x) in (A.3) one can undo this transformation and replace F{x) by the two loop topology of 
figure 7. As the final x-integration is then the same for both sides one can relate the values of the 
integrals on both sides quite straightforwardly. Indeed, expanding out the Gram determinant 
and expressing all the scalar products in terms of factors which appear in the denominator, one 
obtains the general result which is implicit in p2|. We have 



(ai, 02, as, "4, "5) 



[-1,0,0,0,-1] - [-1,-1,0,0,0] + [-1,0,-1,0,0] 

+ [0,0,0,-1,-1] + [0,-1,0,-1,0] - [0,0,-1,-1,0] 

+ [0,-1,0,0,-1] + [0,0,-1,0,-1] - [0,0,0,0,-2] 

- [0,0,0,0,-l]_ - [-1,0,0,-1,-1]+ + [-1,-1,0,-1,0]+ 
+ [-1,0,-1,-1,0]+ + [-1,0,-1,0,-1]+ + [-1,-1,-1,0,0] 

- [-1,0,-2,0,0]+ - [-2,0,-1,0,0]+ + [0,-1,0,-1,-1]+ 

+ [0,-1,-1,-1,0]+ - [0,-2,0,-1,0]+ - [0,-1,0,-2,0]+ 

1 



[0,-l,-l,0,-l]_ 



2(//-l)(2^-l) 



(A.4) 



where 



[ni,n2,n3,n5,n5]„g = (ai + ni, 02 + "-2, "3 + "-3, "4 + ?^4, "5 + "-5) 



(A.5) 



and the left hand side of ( |A.4| ) represents the value of the graph of figure 7 but evaluated in 
(d + 2)-dimensions. 
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